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This work describes a class of topological phonon modes, that is, mechanical vibrations localized at 
the edges of special structures that are robust against the deformations of the structures. A class of 
topological phonons was recently found in 2-dimensional structures similar to that of Microtubules. 
The present work introduces another class of topological phonons, this time occurring in quasi 
one-dimensional filamentary structures with inversion symmetry. The phenomenon is exemplified 
using a structure inspired from that of act in Microfilaments, present in most live cells. The system 
discussed here is probably the simplest structure that supports topological phonon modes, a fact 
that allows detailed analysis in both time and frequency domains. We advance the hypothesis that 
the topological phonon modes are ubiquitous in the biological world and that living organisms make 
use of them during various processes. 

PACS numbers: 63.20.D-,87.10.Hk,87.16.Ka 



Condensed matter research has been profoundly 
marked by the discovery of a new phase of matter 
called the topological insulating phase, 1 4 which has al- 
ready generated a plethora of amazing applications and 
changed the way we understand condensed matter. 5 ' 6 
The hallmark of a topological insulator is the emergence 
of electronic states along any edge cut into a sample of 
such material. Although insulators when probed deep 
inside the bulk, these materials conduct electricity along 
their edges without any resistance. This property can- 
not be destroyed by any chemical or mechanical treat- 
ment of the edge. 7 ' 8 The phenomenon is not restricted 
only to electronic systems. Electromagnetic waves can 
display similar effects and it was shown that specially 
designed photonic crystals can support similarly robust 
electromagnetic wave- modes along their edges. 9 ' 10 Like- 
wise, specially designed mechanical lattices can support 
topological phonon modes. 11 

Edge and surface phonon modes are not rare occur- 
rences in hard matter systems. However, most of these 
modes are very sensitive to the deformations of the sys- 
tems and they can be easily suppressed by various treat- 
ments of the edges or surfaces. There is a special ef- 
fect that occurs in the presence of a magnetic field, the 
phonon Hall effect, 12 ' 13 where robust vibrational edge 
modes can be observed very much like the edge elec- 
tronic states in the Quantum Hall Effect. The effect 
was recently proposed to have a topological nature, 14 by 
following an analysis similar to that of Ref. 11. As op- 
posed to the phonon Hall effect, the emergence of the 
topological phonon modes described in the present work 
(and in Ref. 11) does not require any external field and 
is a consequence of the special intrinsic properties of the 
structures. 

The lattice presented in Ref. 11, exhibiting topolog- 
ical phonon modes, was inspired from the structure of 
Microtubules, 15 a biomaterial synthesized by all living 
cells. Microtubules display a phenomenon known as dy- 
namic instability (DI), 16 19 in which they randomly grow 
and shrink in length, a process that is essential for their 



normal functioning. Most of today's chemo-therapies 
against cancer target Microtubules in an attempt to in- 
hibit their DI. 20 For this reason, understanding the mech- 
anism of DI is one of the most active research areas in 
biomedical research. Ref. 11 showed that the phonon 
spectrum of the dimer lattice of Microtubules displays 
Dirac cones like electrons do in graphene. It also showed 
that these Dirac cones can be split, leading to topologi- 
cal phonon modes, which were indeed observed in explicit 
calculations for ribbon geometries. It was then hypothe- 
sized that the topological edge modes play a major role 
in the DI. 

Topological phonon modes may also be relevant to the 
dynamics of Microfilaments, which control cell motility 
by applying force on the cell membrane. 18 ' 21 They are 
made of the protein actin, arranged in a double-helical 
formation. 22 A pool of ATP-actin monomers is present 
in the cell, and actin polymerization draws on this pool. 
The ATP-monomers collide and bond with the ends of 
the existing filament branches, elongating them. The 
bound ATP-actin hydrolyzes into ADP-actin almost in- 
stantly, releasing quanta of about 12 kT energy. 23 At the 
same time, ADP-actin at the other end of the Microfila- 
ment network depolymerizes and returns to the pool as 
ATP-actin. Eventually capping proteins stop the growth 
of the branches. 21 

One un-explained issue is the way in which the Micro- 
filaments continue to grow while their ends push against 
the cell membrane, as is the case when they generate 
the motile force. 21 It is difficult to explain how the ATP- 
monomers from the solution are still able to collide with 
the ends of the filaments. According to the Elastic Brow- 
nian ratchet model, 21 ' 24 the Microfilaments vibrate as 
spring-like wires and the edges adjacent to the cell mem- 
brane bend laterally, exposing the ends to the pool of 
ATP-actin. This allows additional actin monomers to 
squeeze in and attach themselves to the ends of branches, 
as illustrated in Fig. 1. The restoring force straightens 
the Microfilament, which pushes against the cell wall gen- 
erating the motile force. 
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FIG. 1. (Color online) The ends of the Microfilaments oscil- 
late while pushing against the membrane, exposing their ends 
to the pool of available ATP-Actin. The ATP-Actin is then 
able to attach to the exposed ends and continue the filament's 
polymerization. 



We hypothesize that the bending of the Microfilaments 
is caused by topological edge modes, powered by the 12 
kTs released during hydrolysis of the ATP-actin. The 
present work demonstrates the existence of such modes 
in filamentary structures similar to that of the Microfila- 
ments. As we shall see through explicit simulations, such 
edge modes do not allow the energy to dissipate into the 
bulk of the filaments, and could indeed lead to vigorous 
shakeup of the ends of the structures, even when excited 
with weak stimuli. The phonon modes discussed in this 
work are of a different type from those previously found 
in Microtubules, which required a 2D structure and spe- 
cial interactions. We do not exclude that newly found 
modes may exist in Microtubules. 

The paper gives a general classification of filamentary 
structures with inversion symmetry and provides a simple 
and direct criterium to identify those structures support- 
ing topological phonon modes. The paper also presents 
an explicit mechanical lattice that displays this new phe- 
nomenon, which is investigated through analytic calcu- 
lations and computer simulations in both time and fre- 
quency domains. 



I. A Z 2 TOPOLOGICAL CLASSIFICATION OF 
THE STRUCTURES WITH INVERSION 
SYMMETRY 

The topological classification of the electronic sys- 
tems with inversion symmetry was recently discussed in 
Refs. 25 and 26. Here we will adopt the 1 dimensional 
discussion presented in Ref. 25 to the context of mechani- 
cal waves. Consider a periodic ID harmonic lattice, with 
K degrees of freedom £ per repeating cell, governed by 
the following equations of motion: 

where t's are coefficients specific to each structure (see 
for example Eqs. 23 and 24). The equations of motion 
for any periodic system slightly perturbed from its equi- 
librium configuration can be cast into the form shown in 



Eq. 1. The ansatz of traveling waves: 

C(t) = Re^K^")], k e [-tt.tt] (2) 
leads to the usual normal modes equation: 

[M(k)-uj 2 }A(k) = 0, (3) 

where M(k)a,p = E m C/3 e_i/cm and A k ) is the K di- 
mensional vector encoding the amplitudes A a (k). Let 
oj s (k) and A 3 (k), s= 1, . . . , K, be the solutions to Eq. 3. 
The solutions at fc=±7r are identical. Also, A s (k) are 
assumed to be normalized to one: A*(k) • A 3 (k)=l, for 
all s= 1,...,K. 

One essential requirement is the existence of a gap in 
the vibrational spectrum of the system, that is, an inter- 
val [cj_,cj + ] empty of normal frequencies: u s (k) < uo- 
for s = 1, . . . , S and u) s (k) > cj + for s = S + 1, . . . , K, 
with uj+ strictly larger than uj- . We refer to the interval 
[cj_,cj + ] as the spectral gap. 

For each the K normal modes A s (k) generate a K- 
dimensional complex space C K . The S modes with fre- 
quencies below the spectral gap generate a ^-dimensional 
subspace <S, which we should view as an S'-dimensional 
hyperplane of the if -dimensional complex space. When 
the ^-number is changed, this hyperplane twists, much 
as a Moebius band will in real 3-dimensional space. Our 
task is to classify the twists of the ^-dimensional hyper- 
plane. 

A. General Arguments 

We assume that the ID harmonic lattice has inversion 
symmetry, which implies the existence of a k- independent 
matrix V such that 

VM(k)V- x = M(-fc), (4) 

where V is an unitary matrix that squares to indentity: 

VV ] = I, VV = I. (5) 

The inversion symmetry and the matrix V will be dis- 
cussed in great detail for our explicit example. In the 
following, we demonstrate that the harmonic lattices 
with inversion symmetry fall into at least two topolog- 
ical classes. 

The following analysis is standard in Quantum 
Mechanics, 27 but is definitely not common for the present 
context. We need to define a parallel transport, that is, 
a rule that tells us how to change a vector in the hy- 
perplane S when k is modified so that it stays parallel. 
The construction of the parallel transport is not unique, 
however, the emerging topological classification is inde- 
pendent of how the parallel transport is defined. 28 Let 
P{k) denote the projector onto the S'-dimensional hyper- 
plane generated by the normal modes with frequencies 
below the spectral gap: 

[P(k)U = T, 8 .=i[Mk)UMme- ( 6 ) 
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FIG. 2. (Color online) The adiabatic transport is carried over 
the paths 7 and 7' . 

This is a KxK matrix and ^-dimensional hyperplane S 
is then simply given by P(k)C K . A parallel transport 
can be defined by the so called monodromy U(k, ko), de- 
fined as the unique solution to the equation: 

d k U(k, ko) = i[P{k), d k P(k)]U(k, fe ), (?) 

with the initial condition U(ko,ko)=P(ko). Once we 
computed this U(k,ko) for each ke[— 7r,7r], we can de- 
fine the parallel transport as the map that takes an ini- 
tial vector A(ko)=Yls=i c sA s (ko) of the hyperplane at ko 
into the vector A(k)=U (fc, ko)A(ko) of the hyperplane at 
k. Physically, this parallel transport gives the change in 
a vibrational state of the system when the frequency is 
adiabatically changed. 

Since the normal modes at /c=±7r are identical, we can 
identify these two k points and think that k is defined on 
a circle as in Fig. 2. If we do so, then we see that U (tt, —tt) 
takes the hyperplane P(—tt)C k into himself, becoming 
an unitary matrix that we will call U 11 to relate it to 
the paths shown in Fig. 2. It is known from the non- 
Euclidean geometry that the parallel transport along a 
closed path will not necessarily return a vector into itself. 
Something similar happens here, that is, the map U 1 does 
not return a vector into itself. 

The monodromy U 1 has special properties when the 
inversion symmetry is present. Indeed, assuming ko=— 7r, 
a conjugation of Eq. 7 with V gives: 

dkiPUik,-*)?- 1 } 

(8) 

= i[P(-k), d k P(-k)]PU(k, -n)V~\ 

with the initial condition P/7(— 7r, — 7r)P _1 =P(7r). This 
is just the equation for U{— k,ir), which simply shows 
that VU(k, — 7r)P _1 coincides with U(—k,7r). Equiva- 
lently, we can think that V sends the path 7 into the 
path 7' of Fig. 2. Now obviously U^Uy equals identity, 
therefore: 

detiU^VU^V- 1 } = 1. (9) 

Using the elementary properties of the determinant, 

det{AB} = det{A}det{P}, det{P -1 } = det{P}~\ 

(10) 



we conclude that det{£/" 7 } 2 =l. Consequently, the deter- 
minant can take only two values: 

det{£> 7 } = ±1. (11) 

This is one of our main conclusions. It shows that the 
set of matrices M(k) satisfying Eq. 4 can be split into two 
categories C±, the criterium being the value of det{£7 7 } 
for the corresponding monodromy. A matrix M(k) that 
was placed in one category can not be morphed into a 
matrix from the second category by a continuous defor- 
mation that keeps the spectral gap open. Indeed, under 
such deformation, U 1 will change smoothly and, conse- 
quently, its determinant must change smoothly. Hence, 
it cannot make sudden jumps between +1 and —1. If 
the spectral gap closes, than P(k) becomes ill defined 
and the monodromy U 1 can no longer be defined. If the 
gap opens again, the monodromy U 1 can emerge with a 
different determinant. 

In the following, we describe how one can determine 
the signature of det{/7 7 } through an elementary calcula- 
tion. We have successively: 

det{£> 7 } = det{£7(7r, 0)17(0, -tt)} 
= det{0"(7r, 0)Pt)(0, tt)^ 1 }. ^ 

Taking into account that V~ 1 =V and inserting the ap- 
propriate projectors, we obtain: 

det{£7 7 } = det{P(7r)^(7r,0)P(0)PP(0)[>(0,7r)P(7r)P}. 

(13) 

Using again the elementary properties of the determinant 
and the fact that U(tt, 0)C/(0, 7r)=Id we obtain: 

det{[> 7 } = det{P(0)PP(0)} det{P(7r)PP(7r)}. (14) 

In other words, all we need is to compute the action of in- 
version symmetry operation V at the special points k=0 
and k=7r. As we shall see, for concrete examples this can 
be accomplish through straightforward calculations. 

B. Existence of the edge modes 

We show here that the systems with det{U 7 }=— 1 are 
likely to display edge phonon modes. For this, we imag- 
ine the following experiment. We take the infinitely long 
chain and we synchronously weaken the interactions be- 
tween the 0th and 1st cell, and between the Nth and 
(N + l)th cell, and so on, such that when these inter- 
actions are set to zero, we obtain decoupled finite chain 
pieces of length TV. We will call the action we just de- 
scribed the decoupling process. During decoupling, we 
assume that the inversion symmetry is preserved. 

At any moment of the decoupling process, we can see 
the altered chain as a new periodic chain whose repeating 
cell contains N cells of the original chain. Therefore the 
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FIG. 3. (Color online) The structure consists of a periodic 
array of dimers, connected by four different springs. The di- 
agram depicts the system at equilibrium. 

whole discussion of the last subsection still applies. If one 
computes the monodromy matrix U 1 and its determinant 
for the new chain, he will find the following. At the be- 
ginning of the decoupling process, with the interactions 
untouched, det{Z7 7 }=±l for chains belonging to C± cate- 
gories, respectively. At the end of the decoupling process, 
M(k) becomes independent of k and consequently either 
/7 7 =P(— 7r) and det{C/ 7 }=l, or U 1 cannot be defined be- 
cause the system became gapless. For a system in the 
C- category, both scenarios imply that the spectral gap 
closes during the decoupling. Thus, at least one pair of 
phonon bands must emerge in the spectral gap, one band 
coming from the upper and and one from the lower edge 
of the spectral gap. These two bands must move towards 
each other until they touch. In very special instances, 
which we were actually never able to observe, these two 
bands may continue to move until they disappear from 
the spectral gap. But in general, the bands end up inside 
the spectral gap at the end of the decoupling process, on 
top of each other. The bands become dispersionless, that 
is, completely flat, and the modes corresponding to these 
flat bands represent vibrational modes localized at the 
two ends of each finite piece of chain. This is in striking 
contrast with the case of simple harmonic lattices, that 
is, periodic arrays of harmonically interacting (point like) 
bodies without internal structure, where it was concluded 
in Ref. 29 that vibrational modes localized at the edges 
are very unlikely to occur in one dimensional structures. 

The vibrational modes discussed above are robust, in 
the sense that, no matter how the decoupling is done, 
these edge modes will always emerge during the decou- 
pling process. For example, a real lattice might become 
strongly distorted near the edges when a finite piece is 
being cut out but, even in this case, we still expect to see 
edge phonon modes for a system in the C_ category. 



II. EXAMPLE OF A TOPOLOGICAL STRING 

In this section we introduce our example. It is the 
structure shown in Fig. 3, made of an array of dimers, 
consisting of two masses M linked by a rigid and massless 
arm of length 2d. The adjacent dimers are connected 
by four different ideal springs. We will allow only two 
degrees of freedom per dimer: a horizontal motion of 
the center of mass, whose position will be marked by x n 
(y n will be constrained to 0), and a pivotal motion of the 



FIG. 4. (Color online) If K++=K--, the structure is mapped 
into itself by the symmetry operations shown in this diagram. 

dimer in the plane of the structure and around the center 
of mass, described by the angle <j) n . We assign the labels 
± to the masses on the upper/lower rows of the structure. 
Using these labels, the springs between two dimers can 
be uniquely labeled by a/3, depending on which masses 
are connected by it. The corresponding spring constants 
will be denoted by K a p. 

Let us point out that if we chose K ++ =K , then 

the system is symmetric to the operations described in 
Fig. 4. These operations define the inversion symmetry 
of our system. 

A. The Equations of Motion 

We decided to discuss the equations of motion for the 
system in detail, in order to assist those who want to ex- 
plore beyond what is presented here. With the notations 
introduced in Fig. 5, the positions of two masses of the 
n-th dimer can be conveniently written as: 

= x n — adsin((j) n )] y% = adcos((j) n ). (15) 

The equations of motion for the n-th dimer are: 

2Mx n = F n , x , 2M^ n = r n , (16) 

where F UjX and r n represent the net horizontal force and 
the net torque on the n-th dimer. They are given by: 

and 

T n = S [( x n ~ X n)(Fn,y ~ ^n-l,y) 

<*,P (18) 

cfnK n,x n—l,xJ] 

where F^f x , are the horizontal/ vertical components of 
the forces generated by the springs at the right side of 
the n-th dimer: 

F?P = K aP (l - lf/lf)(y? n+1 - ft), 

with: 

If = V^n+l-^ + O/f+l"^) 2 (20) 

being the distance between the a and f3 masses of the 
n-th and n + 1-th dimers, and 1^ being the un-stretched 
length of the af3 spring. 
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FIG. 5. (Color online) The diagram introduces the notations 
used in the text. The array of dimers is shown at the equi- 
librium (faded) and at some arbitrary snapshot in time. x n 
represents the position of the center of mass of the n-th dimer 
and 4> n represents the angle of rotation of the n-th dimer rel- 
ative to the vertical axis. In the equilibrium configuration, all 
(j) n are equal to <po. 
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form of the equations of motion: 

= — | S ^afii^n+l — 2£ n + £n-l) 

a(3 



(23) 



\ E cos O + K V p sin 0o] + 

a/3 



and 

-'(I = - E cos O + «S» sin 0o]^ 

Q! 

+ 5 E "[O^ cos 0o + 0^ sin O ] + 
+ E { cos2 + fi^ sin 2 O ] + 

~ E cos 2 0o + SQ sin 2 0o](C 2 +1 + e 2 _i)- 

a/3 

(24) 

The coefficients have the following explicit expressions: 



(25) 



FIG. 6. The phonon bands for (a) K = K++ = 0.5 and 

K-+ = K+- = 2 and (b) K— = K ++ = 0.5 and K-+ = 
K+- = 4. 



The above equations of motion have been implemented 
numerically and the time evolution of the structure has 
been investigated using a 4th order Runge-Kutta time 
propagator. The results will be presented shortly. In all 
our calculations, we used length and time units so that 



M = 1 



In 



1 and d was set to 1. 



B. The Small Oscillations 

To explain and understand the features seen in our 
time domain simulations, it is useful to compute first the 
phonon spectrum of our structure, that is, the relation 
between the frequency and the wavenumber of the wave 
propagating modes. Let us introduce the notation: 



t 1 



• r eq £ 



(21) 



In the linear regime, that is, when the dimers are only 
slightly perturbed from their equilibrium configurations, 
we have: 



■ + 4n - OL COs((/> )£n> Vn = ~ a COs(0 O )fn- (22) 



Plugging these expressions in Eq. 16 and retaining only 
the linear terms in ^ and £ 2 , after tedious but otherwise 
straightforward calculations, we obtained the linearized 



{' 



( &V a p x 2 



} 



and 



otp ap 



(26) 



(27) 



with Ax a p=x%+ 1 -xP and Ay a/3 =y% +1 -yP. The dimers 
are assumed in their equilibrium configuration in these 
last two equations. 

To simplify the equations, we assume from now on that 
none of the springs are under tension at the equilibrium 
configuration and that the 1^ are chosen such that 4>o=0. 
The ansatz Cn= Re { A i ei(ut ~ kn) } leads to the following 
normal modes equation M(k)A = oj 2 A, with: 

-01 (1 - COS k) J2 a aK aa + °2 sin k £) Q aK a-a ( 28 ) 

-as cos k^2 a aK^_ a , 
where oY s are the Pauli's matrices: 

»i = (? J) ,** = (? "o ) ,»i = (J • (29) 

The normal modes equation has two solutions ui^(k) 
and Ai^(k). 

Let us start the discussion of these solutions from the 

special case when K+j r =K and =K |_. In this 

case, the normal modes equation becomes: 



(30) 
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FIG. 7. The spectral gap opening for (a) K = K++ = 0.5, 

K-+ = 3 and K+- = 5 and (b) K— = 0.3, K ++ = 0.7 and 
K-+ = K + - = 4. 

with 

Cl =Ea i 9 n S[ 1 " C0S *]' (31) 
C2=Ea i 9 n S?[ 1 " a ^ COSfc ]- 

The phonon modes decouple into a mode (corresponding 
to ei) that involves displacements of the centers of mass 
but no pivotal motion and a mode (corresponding to 62) 
which involves pivotal motion and no displacements of 
the centers of mass. 

Depending on the common values given to and 
K__ and to and K |_, the phonon bands can as- 
sume two distinct behaviors as shown in Fig. 6. We are 
primarily interested in the case (b) because, as exempli- 
fied in Fig 7, when the bands cross each other in that 
way, a small mismatch between and K |_ or be- 
tween and K-- will open a gap in the spectrum, 
thus realizing the last requirement of our general theory. 

C. The Edge Spectrum 

As already mentioned, if we set K++=K , the sys- 
tem becomes symmetric to the inversion operation shown 
in Fig. 4 and the M(k) matrix becomes: 

M (*0 = [ Ea(3 K% - E a COS k] I 

+fr 2 Sin k E a aK a-a - ^3 COS k E a aK a-a' 

In the /c-space, the inversion symmetry is implemented 
by V = cr 3 and one can explicitly verify that: 

a 3 M(k)a 3 = M(-fe). (33) 

Moreover, the projectors at the special points k = and 
7r can be computed explicitly. Assuming ^ a aK%- a < 0, 
they are: 

P(0) = \{I - a 3 ), P(tt) = |(7 + a 3 ) (34) 

If J2 a a ^a-a > 0? ^ ne expressions will be switched. In 
both cases, we can explicitly verify that: 

det{P(0)PP(0)} det{P(7r)PP(7r)} = -1, (35) 



and consequently the system is topologically nontrivial 
and we should observe edge modes. Note that when com- 
puting the determinants in Eq. 35, one needs to discard 
the null eigenvalues. 

If we keep K |_ = and consider a difference be- 
tween and the system still has an inversion 
symmetry but in this case the symmetry is implemented 
in the /c-space by the identity matrix. Hence: 

det{P(0)PP(0)} det{P(7r)PP(7r)} = +1. (36) 

Consequently, the system is trivial and we should not 
observe robust edge modes. 

To verify these predictions, we performed the follow- 
ing numerical experiment. We considered a chain of 100 
dimers linked in a circle configuration (the radius of the 
circle is very large so any physical bending of the struc- 
ture can be ignored) . In our theoretical argument of Sec- 
tion IB we used a gedanken experiment involving an in- 
finite chain which was gradually sectioned in finite equal 
pieces. In practice, we cannot simulate such an infinite 
chain, but we can still start from a configuration with no 
edges (to avoid un- wanted phonon reflections), which is 
precisely the circle configuration mentioned above. Then 
we slowly weakened all four spring constants between two 
dimers, by making the substitution K — ^ wK with w con- 
tinuously varying from 1 to 0. When w = 0, the circle is 
fully opened in a finite piece with two separated edges, 
exactly like in the gedanken experiment. If topological 
phonon modes are present, we expect to see them grad- 
ually emerging from the bulk spectrum as described in 
Section IB. 

For each w we have computed all 200 normal modes 
and we placed their frequencies on a vertical axes. Fig. 8 
shows the results of these calculations, which illustrate 
how the frequencies of the normal modes change as w 
is varied from 1 to 0. In the topological case = 
K++ — 0.5 and = 3 and K |_ = 5 one can ob- 
serve two solitary normal frequencies separating from the 
bulk spectrum and moving towards each other. When 
the two edges are completely formed at w = 0, the two 
frequencies meet near the middle of the gap (the oscilla- 
tory motion of these modes is discussed in the next sec- 
tion). In contradistinction, no such modes are observed 

for the non-topological case K = 0.7, = 0.3 and 

K |_ = = 0.5. This explicit calculation confirms 

our general theoretical predictions. 



D. Time Domain Analysis 

As already mentioned, the (non-linear) equations of 
motion have been implemented numerically and the mo- 
tion of the dimmers has been studied in real time 30 . Here, 
we discuss the manifestation of the edge modes in the real 
time dynamics of the dimer chain. For this, we consid- 
ered a chain containing 7V=100 dimers plus two addi- 
tional "edge" dimers. We fixed the very left dimer to the 
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FIG. 8. Each vertical sequence of dots represent the frequen- 
cies of the normal modes of a chain made of 100 dimers ar- 
ranged in a circle. The spring constants K of the springs con- 
necting two dimers were gradually weakened by making the 
substitution K^wK and letting w vary from 1 to 0. The fre- 
quencies were recomputed for many w values. Panel (a) refers 
to the topological case, and panel (b) to the non-topological 
case. 



upward position while forcing the very right dimer into 
the following motion: 



^at+i = Ai shiest, (/>at+i = A2 sincjt. 



(37) 



The rest of the N dimers were released from the equilib- 
rium with zero velocities (see Fig. 9). 

The resulting motion of the dimers has been plotted 
every At =3, and movies were generated for various pul- 
sations uo. For a more complete representation, we chose 
to look at (1) the actual dimer positions, (2) only at the 
horizontal displacement x n (t) of each center of mass, and 
(3) only the angular displacement <p n (t) of the dimers as 
a functions of time. In the following, we will make refer- 
ence to the wave modes discussed in subsection IIC. The 
movies (see the supplemented material) reveal that, in- 
deed, if uo is below approximately 1.4 (see Fig. 6), only the 
wave propagating mode that involves the translation of 
the dimers is excited, even though A2 was given the same 
value of 0.03 as A\. Above 1.4 and below the edge of the 
spectral gap we see both wave modes being excited. Sim- 
ilar behaviors are observed above the spectral gap. What 
really interests us is what happens when uo takes values 
inside the spectral gap. In order to generate a mean- 
ingful plot, we let the movies progress without erasing 
the images already played out. The resulting plots shows 
the trajectories followed by each dimer and, in particu- 
lar, they reveal the amplitudes of the oscillating motion 
for each dimer at the time when the movie was stopped. 
In Fig. 9, the last dimer was forced into the motion de- 
scribed in Eq. 37 with Ai=A2=0.3 and co>=1.6, chosen 
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FIG. 9. (Color online) Top diagram: Representation of the 
dimer chain and the forced motion of the rightmost dimer. 
Main panel: The trajectories covered by the dimers after 
t=500, 1000 and 2000. Each sub-panel contains three plots, 
representing, from the top to bottom, the actual trajectories 
of the dimers, the trajectory of the horizontal displacements 
x n (t), and of the angular displacements <f> n (t). Plots (a) (b) 
and (c) represent the behavior of the dimers in the topological 

case where K — K-\--\- — 0.5 and K-\ — 3 and K |_ = 5. 

Plots (a') (b ; ) and (c) represent the behavior of the dimers 
in the non-topological case where K-- = 0.7, = 0.3 and 

K-+ = K+- = 4. 



to be in the middle of the spectral gap. The movie was 
allowed to run for three time intervals: 500, 1000, 2000. 

The left panels in Fig. 9 simulate the situation illus- 
trated in Fig. 8(a). This is the topological case, for which 

the spring constants were set to K = = 0.5 and 

K + _ — 3 and K h = 5 and an edge resonance is ex- 
pected. As one can see, only the amplitudes of the dimers 
close to the right edge reach appreciable values and the 
amplitudes are seen to decay exponentially away from 
the edge, demonstrating that we are indeed dealing with 
an edge mode. Furthermore, one can see the amplitudes 
increasing as time progressed, the angular displacement 
of the second right dimer reaching an amplitude of ap- 
proximately 0.4 after t = 2000, more than 10 times the 
amplitude of the forced oscillation imposed on the first 
dimer. This reveals an important property of the edge 
resonance, namely the ability to absorb and store energy 
in the proximity of the edge of the structure. 

The right panels in Fig. 9 simulate the non-topological 
case presented in Fig. 8(b), where the spring constants 
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were set to K__ = 0.7, K++ = 0.3 and K_+ K+_ 
1.0. In this case we don't expect an edge mode and in- 
deed the amplitudes of the dimers near the edge remain 
practically zero at all times. 

There is a small artifact in the plots of Fig. 9. Due 
to a non-stationary effect steaming from the relatively 
large amplitude used for shaking the end dimer, the av- 
erage positions of the centers of the dimers slowly drift to 
the left in time. Because of the way we generate Fig. 9, 
this drift may give the impression of a finite oscillation 
amplitude for dimers away from the edges. We prompt 
the reader to watch the movie, where one can clearly see 
that the oscillation amplitudes of the dimers away from 
the edges are practically zero (the slow drifting motion 
is also visible). 

III. CONCLUSIONS 

Using a monodromy argument, we have demonstrated 
that the filamentary structures with inversion symme- 
try fall into at least two topologically distinct classes. 
We have argued that one of the two classes contain sys- 
tems that should display robust topological edge modes. 
The general theoretical predictions were verified using an 
explicit example of a mechanical structure that display 
robust edge modes. 

The structure was inspired from that of actin Microfil- 



aments. We hypothesized that the energy from the ATP 
hydrolysis in Microfilaments is stored in these robust vi- 
brations and is then used to drive the motion described 
in the Elastic Brownian ratchet model. This opens an 
interesting research direction, which is worthwhile pursu- 
ing because we believe that, with the present technology, 
well designed and focused experiments can reveal if the 
Microfilaments have or have not edge modes. 

Although not very accurate, the model allowed us to 
present the concept of topological phonon modes in a 
very explicit and detailed exposition, without being de- 
railed by unnecessary complications. The simplicity of 
the structure will allow detailed analyses of other inter- 
esting questions, such as what happens when defects are 
present throughout the bulk of the chain, how is the edge 
excited by the hydrolysis of GTP actin or tubulin, how 
is damping effecting the picture, etc.. The structure pre- 
sented here can become a very useful pedagogical tool to 
introduce the concept of topological phonon modes in an 
accessible and explicit way via computer simulations or 
real lab observations. 
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